library(tidyverse)
library(sf)
library(tigris)
library(censusapi)
library(leaflet)

Sys.setenv(CENSUS_KEY="b88495f8315f45b2d100dee9ba8f4c489a7371c2")
blockgroups <- block_groups("CA",037,progress_bar = F)

# list <- c("060372060102","060371997001","060371997002","060379800101","060371977003","060372071032","060372060101")
# 
# narrow_blockgroups <-
#   blockgroups %>% 
#   filter(GEOID %in% list)
  
 test <- getCensus(
      name = "acs/acs5",
      vintage = 2018,
      region = "block group:*",
      regionin = paste0("state:06+county:037"),
      vars = "B01003_001E"
    ) %>%
  transmute(
    origin_census_block_group = paste0(
      state,county,tract,block_group
    ),
    pop = B01003_001E
  )


test1 <- test %>% left_join(blockgroups,by = c("origin_census_block_group" = "GEOID")) %>% dplyr::select(origin_census_block_group, pop,geometry) %>% st_as_sf()
leaflet() %>% 
  addTiles() %>% 
  addPolygons(
    data = test1,
    fillColor = "blue",
    color = "white",
    opacity = 0.5,
    fillOpacity = 0.5,
    weight = 1,
    label = ~pop,
    highlightOptions = highlightOptions(
      weight = 2,
      opacity = 1
    )
  )